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Abstract 



The possible discovery of the + pentaquark has motivated a number of studies of its nature 

O: 

using lattice QCD. While all the analyses thus far have focused on spin-^ states, here we report 
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the results of the first exploratory study in quenched lattice QCD of pentaquarks with spin i. For 



the spin- 1 interpolating field we use a product of the standard N and K* operators. We do not 



find any evidence for the standard lattice resonance signature of attraction (i.e., binding at quark 

masses near the physical regime) in the J p = | channel. Some evidence of binding is inferred 

in the isoscalar | channel at several quark masses, in accord with the standard lattice resonance 

signature. This suggests that this is a good candidate for the further study of pentaquarks on the 

Qh. lattice. 
43 
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I. INTRODUCTION 



The recent reported observations of the strangeness +1 pentaquark, G + , having minimal 
quark content uudds, have led to a tremendous effort aimed at understanding its properties 
both experimentally and theoretically. Many model studies and phenomenological analyses 
have explored various aspects of its structure and production mechanisms, and have at the 
same time have revealed a number of challenges for its interpretation as an exotic resonance 



DDDD 



with a particularly narrow width (for recent experimental reviews see Refs 

To tackle this problem from first principles in QCD, a number of lattice studies have 
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12. U 



14 



15. 161. These have used various 



recently been undertaken 

local spin-| interpolating fields (either A^-type or diquark-diquark-s type), and, in the 
15| . also non-local fields. Several of these studies have interpreted their results 



case of Ref. 



as indicating the presence of a resonance, while others report signals which are consistent 
with NK scattering states. 

A major challenge in the lattice studies has been the identification of a resonance state 



from the NK scattering states. Several groups have sought to distin. 



and scattering states by comparing the masses at different volumes 



mish 



3, HQ QQ 



;ne resonance 
The 



volume dependence of the residue of the lowest lying state has also been proposed as a way 
to identify the nature of the state ^, w| . Alternatively, hybrid boundary conditions have 



been used in Refs. 



101 ] to differentiate the resonance in the negative parity channel from 



the S'-wave NK scattering state. 

In Re, Q we e m pl„ y ed a co m pli m en tary approach to invoke S pi„-| penta^ 
resonances by searching for evidence of sufficient attraction between the constituents of 
the pentaquark state such that the resonance mass becomes lower than the sum of the free 
decay channel masses. We labeled this pattern as "the standard lattice resonance signature" 
because this signature is observed for conventional baryon resonances studied on the lattice 

ay y y&. By m ^ _ of ^ « ^ t0 the mass of 

the decay channel we found no binding at any quark mass and hence no evidence for such 
attraction. The absence of binding cannot be used to exclude the possibility of a resonance, 
as the attractive forces simply may not be strong enough to provide binding. On the other 
hand, the presence of binding would provide a compelling resonance signature warranting 
further study. 
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One of the major puzzles in pentaquark phenomenology has been the anomalously small 
width (< 1 MeV) observed in the experiments which have produced a positive signal. A 
possible explanation for this may be that if the pentaquark has J p = | , its decay to N + K 
must be via a D-w&ve, which would consequently be suppressed. In this paper we therefore 
extend the analysis of Ref. p] to spin-| pentaquarks. We examine both the positive and 
negative parity states, in both the isoscalar and isovector channels. In Sec. II we describe 
the interpolating field and outline the lattice techniques employed in this analysis. The 
results are presented in Sec. Ill, where we discuss in detail the mass splittings between 
the pentaquark and two-particle scattering states. Finally, conclusions and suggestions for 
future work are summarised in Sec. IV. 



II. LATTICE DETAILS 



A. Interpolating fields 

The simplest iVi^-type interpolating field used in lattice simulations, referred to in Ref. 
as the "colour-singlet" NK field, has the form: 

Xnk = -±=e abc (u Ta C l5 d b ) {u c (s e t l5 d e ) T (u <- d)} , (1) 

where the — and + corresponds to the isospin 1 = and 1 channels, respectively. This field 
has spin |, and transforms negatively under the parity transformation q — > 7 g. 

One can access spin-| states by replacing the spin-0 fT-meson part of Xnk with a spin-1 
K* vector meson operator, 

Xnk* = ^=e abc (u Ta C l5 d b ) {u c {s^d e ) T (u <- d)} , (2) 

where again the — and + corresponds to the isospin 1 = and 1 channels, respectively. The 
field Xnk* transforms as a vector under the parity transformation, and has overlap with 
both spin-| and spin-| pentaquark states. States of definite spin can be projected from 
Xnk* by applying appropriate projectors, as discussed in the next section. 
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B. Lattice Techniques 



The masses of the spin-| and spin-| pentaquark states are obtained from the two-point 
correlation function 

Gnv{t>P) = ^2 exp(-ip-x) {0\T Xn{x) Xu(0) I 0) . (3) 

X 

To project states of definite spin from the correlation function Q^ v it,'p) we apply the spin 
projection operators |3| 

1 1 

PUV) = 9,u ~ pUp) , (4) 

for spin- 1 and |, respectively. 

The spin-projected correlation function receives contributions from both positive and 
negative parity states. The use of fixed boundary conditions in the time direction enables 



states of definite parity to be projected using the matrix 



18, 



22 



r, = If 1± ^ 7l| , ( 5 ) 



2 V E B ± 

for negative and positive parities, respectively. We note that this differs from that of 
Ref. [19], where interpolating fields transforming as pseudovectors, in accord with the Rarita- 
Schwinger spinor-vectors, were used. Masses of states with definite spin and parity can then 
be obtained from the spinor trace of the spin and parity projected correlation functions, 

G(t,p) = tT sp [rg, u (t,p)P^( P )] 

= X B X B exp (—Est) 

B 

*=°° A A°exp(-m t) . (6) 

This function is a sum over all states, B, with energy Eg, and X B and X B are the couplings 
of the state B to the interpolating fields at the source and sink, respectively. Note that as 
in Ref. [3], we consider the case /x = 3 to reduce the computational cost of the calculation. 

Since the contributions to the two-point function are exponentially suppressed at a rate 
proportional to the energy of the state, at zero momentum the mass of the lightest state, 
mo, is obtained by fitting a constant to the effective mass, 

G{t, 0) 



M etl (t) = In 



G(t + 1,0) 
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2 ) A(1232) and the P-wave N+tt decay channel. 



t— >oo 



m . (7) 

Following our previous work , we search for evidence that the resonance mass has become 
smaller than the sum of the free decay channel mass for pentaquark states created by the 
interpolating field Xnk*- For this purpose it is useful to define an effective mass splitting. 
For example, in an S"-wave decay channel, 



AM ett (t) = Mg(t)- (Mf(t) + M£(t)) 



m 5q - (m B + m M ) 



(8) 



where Mg ff (t) and M^(t) are the appropriate baryon and meson effective masses for a 
specific channel. For a P-wave decay channel, the effective masses are combined with 
the minimum nontrivial momentum on the lattice, 2tt/L, to create the effective energy, 



E eS (t) = a/ (M cS (t)) 2 + (2n/L) 2 , for each decay particle, where L is the lattice spatial 
extent. The advantage of this technique is that it measures a correlated mass difference, 
thereby suppressing the sensitivity to systematic uncertainties (such as using different fit- 
ting ranges). Moreover, correlations in the effective masses can cancel, leading to a more 
accurate determination of the mass splitting. 

As an example of a lattice signature of a resonance, we consider the mass splitting between 
the J p = 3/2 + A baryon and the energy of the P-wave N + tt decay channel. In Fig. ^ the 
typical resonance signature is clearly seen, where the A is bound on the lattice at heavier than 
physical quark masses, because of its lower energy compared with the free decay channel. 
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TABLE I: Lowest energy decay channels for each pentaquark state on the lattice, where the A 
baryon is bound. 

I(J P ) Decay channel 
0,1(1") 5-wave N + K 
0, l(i + ) P-wave N + K 

0(|") S-waveN + K* 

0(| + ) P-wave N + K 

1(|") S-waveA + K 

1(|+) P-wave N + K 

For a pentaquark resonance we shall apply the same criteria, and consider the mass 
splittings between the pentaquark state and the corresponding baryon and meson free two- 
particle scattering states. The signal we are searching for is evidence of a pentaquark bound 
state at quark masses near the physical regime. In Table we summarise the lowest energy 
decay channels for the various isospin, spin and parity quantum numbers considered in this 
analysis. 

In the case of the A baryon, the binding is seen to become stronger at larger quark 
masses. Indeed, from their minimal quark content in the heavy quark limit, one expects to 
recover a A to N ix mass ratio of 3/5. In the case of a pentaquark resonance, the analogous 
mass ratio will be 1 in the heavy quark limit and the mass splitting will vanish relative to 
the hadronic mass scale in the heavy quark limit. Hence a lattice resonance signature for a 
pentaquark state is binding (a negative mass splitting) at intermediate quark masses, above 
the physical regime, with a general trend of binding as a fraction of hadron mass towards 
zero as the heavy quark limit is approached. 



C. Lattice Simulation Formalism 

This analysis is based on the ensemble of 290, 20 3 x 40 SU (3) gauge-field configurations 
considered in p. Using the mean-field 0( a 2 ) -improved Luscher-Weisz plaquette plus rectan- 
gle action J^] , the gauge configurations are generated via the Cabibbo-Marinari pseudoheat- 
bath algorithm with three diagonal SU(2) subgroups looped over twice. The simulations are 
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performed using a parallel algorithm with appropriate link partitioning, as described in 



24j . The lattice spacing is 0.128(2) fm, determined using the Sommer scale tq = 0.49 



Ref. 
fm. 

For the fermion propagators, we use the FLIC fermion action |25(, an (9(a)-improved 
fermion action with excellent scaling properties providing near continuum results at finite 
lattice spacing j^J. 

A fixed boundary condition in the time direction is implemented by setting Ut(x,N t ) = 
V x in the hopping terms of the fermion action. Periodic boundary conditions are imposed 
in the spatial directions. To explore the effects of the fixed boundary condition we have 
examined the effective mass of the pion correlation function and the associated XcLf obtained 
in various fits. The pion is selected as it has the longest correlation length and will be a 
worst case scenario for the boundary effects. We find that the fixed boundary effects are 
completely negligible prior to time slice 30, which is the limit of signal in the pentaquark 
correlation functions presented below. 

Gauge-invariant Gaussian smearing [2^ in the spatial dimensions is applied at 
the fermion source at t = 8 to increase the overlap of the interpolating op- 
erators with the ground states. Six quark masses are used in the calcula- 
tions, with k = {0.12780,0.12830,0.12885,0.12940,0.12990,0.13025} providing am n = 
{0.540,0.500,0.453,0.400,0.345,0.300} Q The strange quark mass is taken to be the 
third largest (« = 0.12885) quark mass. This k provides a pseudoscalar mass of 697 MeV 



which compares well with the experimental value of a/2M^ — M% = 693 MeV motivated by 
leading order chiral perturbation theory. The error analysis is performed by a second-order, 
single-elimination jackknife, with the \ 2 P er degree of freedom obtained via covariance ma- 
trix fits. Further details of the fermion action and simulation parameters are provided in 
Refs. 3, Q| and [^j respectively. 



III. RESULTS 



In this section we present our results for the masses of spin-| pentaquarks for both 
negative and positive parity, in both the isoscalar and isovector channels. 
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A. Negative parity isoscalar channel 



We begin the discussion of our results with the isoscalar, negative parity channel. The 
effect of the spin projection on the correlation function is highlighted in Fig. |21 This figure 
shows the effective mass plot of the G 33 = tr sp {T~G 33 } component of the correlation function, 
which contains a superposition of spin-- and spin-| contributions. Upon spin projection, as 
described by Eq. (0), two distinct states are identified. 

The effective masses of the spin-| and spin-| states are presented in Figs. El and 
respectively. To extract the masses of the lowest energy states from these effective masses, 
we fit over the time slices 20 — 30 for the spin-| and 18 — 21 for the spin-| states, respectively, 
where these intervals have been selected so as to obtain acceptable values of the covariance- 
matrix based \ 2 P er degree of freedom (Xdof)> which we restrict to x\ f < 1-5- The resulting 
masses are presented in Fig. |3J together with the mass of the I(J P ) — 0(~ ) state extracted 
with the standard NK pentaquark operator 

[5] 

of Eq. and the relevant (non- interacting) 
two-particle states. We obtain the expected result that the masses of the I(J P ) = 0(| ) 
state extracted with the NK and NK* interpolators are in excellent agreement. The mass 
of the I(J P ) = 0(| ) state is also similar to that of the N + K* two-particle state, but 
lies consistently above the latter, suggesting the presence of some repulsion in this channel. 
Therefore we cannot conclude any evidence of a bound state in this channel. 



B. Positive parity isoscalar channel 

Next we consider the isoscalar state in the positive parity channel. Contrary to the 
negative parity signal, the spin projection shown in Fig. |U] has a less pronounced effect on 
the effective masses. The effective masses of the spin-^ and | states are presented in Figs. [7| 
and |Hl respectively. The quality of the signal in the positive parity sector is significantly 
reduced relative to the negative parity channel, as in the spin-| analysis in Ref. 

The effective mass of the spin-| state is fit at time slices 18 — 20 and the spin-| state at 
time slices 19 — 24. The poor quality of the signal limits the analysis to the four largest quark 
masses considered. In Fig. El we show the fitted masses of the two spin states extracted with 
the NK* interpolator. For comparison, we display the mass of the spin-| state extracted 
with the NK interpolator, and the energies of the relevant two-particle states. Once again, 
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FIG. 2: Effective masses of the negative parity projected, isoscalar correlation functions calculated 



with the NK* pentaquark interpolator, Xnk*- ^he mass corresponding to the unprojected G 
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correlation function (squares) is compared with that of the spin-i (triangles) and spin-| (circles) 
projected correlation functions. The data corresponds to our heaviest quark mass. 
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FIG. 3: Effective mass of the I(J P ) = 0(^ ) pentaquark calculated with the NK* pentaquark 
interpolator. The data correspond to ~ 830 MeV (squares), 700 MeV (circles), and 530 MeV 
(triangles). 

we see excellent agreement between the masses of the spin-| states extracted with the NK 
and NK* interpolators. 

Interestingly, the mass of the spin-| state becomes smaller than the non-interacting two- 
particle energy of the P-wave N + K state for intermediate quark masses, i.e., we observe 



1.5 



10 



15 



20 



25 



30 



t 



FIG. 4: As in Fig. El but for the I{J P ) = 0(f) state. 
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FIG. 5: Masses of the I{J P ) = 0(| ) and 0(| ) states extracted with the NK* interpolating 
field as a function of m%. For comparison, we also show the mass of the I(J P ) = 0(i ) state 
extracted from the NK pentaquark interpolator 5], and the masses of the .S-wave N + K and 
N + K* two-particle states. Some of the points have been horizontally offset for clarity. 



binding. As discussed in the previous section and in Ref. [5|, the transition of a resonance 
to a state which lies below the free particle decay channel at quark masses near the physical 
quark masses is the standard resonance signature in lattice QCD. 

Moreover the approach to the heavy quark limit is in accord with expectations. Recall 
that in the case of the + , which has a "fall-apart" decay mechanism, quark counting 
indicates the G + to N + K mass ratio will approach 1 as the heavy quark limit is approached. 



10 



2.0 



1.5 



NK - No Proj. x 

NK* I(J p )=0(3/2+) Proj. 

NK* I(J p )=0(l/2+)Jproj. 
_1 i i i i L i i i i 1_ 



10 



15 



20 



25 



30 



FIG. 6: As in Fig. [21 but for the isoscalar positive parity channel. 
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FIG. 7: Effective mass of the I(J P ) = 0(^ + ) pentaquark obtained from the NK* interpolating 
field. The data correspond to m n ~ 830 MeV (squares), 700 MeV (circles), and530MeV (triangles). 

At intermediate quark masses, however, one expects the resonance signature analogous 
to the A baryon in Fig. ^ and at the two smallest quark masses shown the pentaquark lies 
below the scattering state, which is the necessary condition for the presence of binding. As 
this result presents the possible existence of a pentaquark resonance in the physical quark 
mass regime, it is essential to consider a mass splitting analysis of the effective masses. 

As an indicative example, we review the analysis of the lightest quark mass presented in 
Fig. El In Fig. El we present the effective mass splitting between the five-quark state and 
the lightest non-interacting two-particle state having the same quantum numbers. 

At first sight it might be tempting to consider fits early in the Euclidean time evolution 
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FIG. 8: As in Fig.Q but for the I(J P ) = 0(| + ) state. 
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FIG. 9: Masses of the I(J P ) = 0(^ + ) and 0(| + ) states determined from the NK* interpolating 
field as a function of m^. For comparison, we also show the mass of the I(J P ) = 0(^ + ) state 
extracted with the NK pentaquark interpolator |5|, and the masses of the P-wave N + K and 
N + K* and S-wave N* + K two-particle states. Some of the points have been offset horizontally 
for clarity. 



where the errors are small and a possible plateau catches the eye. However, it is easy 
to demonstrate that such a fit does not describe the lowest lying state in this correlation 
function. For example, Fig. ^2 reports the Xdof f° r a selection of fits to the effective mass 
shown in Fig. where the lower bound of the fit window is fixed at t = 12 (four time 
steps from the source at t — 8), and the upper bound of the fit window is plotted on the 
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FIG. 10: Effective mass splitting between the I(J P ) = 0(| + ) state extracted with the NK* 
interpolator and the energy of the P-wave N + K two-particle state for the lightest quark mass 
shown in Fig. 

horizontal axis. As soon as time slice 19 is included in the fit the Xdof becomes very large 
and continues to increase with the inclusion of time slices 20 and 21. One must therefore 
conclude that the results for time slices 19 through 21 are signal rather than noise, reflecting 
the true ground state of the correlator. The effective mass splitting with a "double plateau" 
as in Fig. El can occur when the interpolating fields couple strongly to a more massive state 
and relatively weakly to the ground state. The former dominates at early Euclidean times 
and the latter at later times. 

To determine an appropriate lower bound for the fit, we begin by returning to Fig. ITU1 
The fluctuation at t = 25 suggests that there is valid data up to time slice 24, after which 
noise begins to hide the signal. Indeed, the Xdof * s approximately invariant for fits of t — 19 
through 24 and beyond. Therefore, time slice 24 is selected for the upper bound of the fit 
window. 

The lower bound of the fit interval must be selected with regard to the systematic time 
dependence of the effective mass and the xiot °f the fit. Figure H~2l reports the latter criteria, 
illustrating the Xdof f° r a selection of fits to the effective mass shown in Fig. El where the 
lower bound of the fit window is shown on the horizontal axis and the upper bound of the 
fit window is fixed at time slice 24. The mass splittings extracted for these fits are shown 
in Fig. El Again, the lower bound of the fit window is shown on the horizontal axis and 
the upper bound of the fit window is fixed at time slice 24. The most statistically precise 
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FIG. 11: The Xdof ^ or a ser i es °f possible fits with a lower bound fixed at time slice 12 and an 
upper bound shown on the horizontal axis. 
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FIG. 12: The Xdof f° r a ser i es °f possible fits with an upper bound fixed at time slice 24 and a 
lower bound shown on the horizontal axis. 

estimate of the mass splitting is obtained for t min = 19, and this result agrees with all other 
determinations at the la level, with the exception of t m j n = 16. 

Having determined that there is genuine signal in the effective mass splitting at time slice 
19, and given the systematic drift in the results approaching time slice 19 as illustrated in 
both Figs. El and EH one has to conclude that nontrivial contributions from more massive 
excited states are still present at time slice 16. The more conservative evaluation is to delay 
the fit to later Euclidean times where the Xdoi ~ 
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FIG. 13: Mass splitting extracted from a series of possible fits with an upper bound fixed at time 
slice 24 and a lower bound shown on the horizontal axis. 

These Xdof analyses have been repeated for all six quark masses considered with similar 
results. Figure H3I shows the mass splitting between the isoscalar spin-| state extracted with 
the NK* interpolator and the P-wave N + K energy, where we fit at time slices 19 — 24 
as concluded above. Here we are able to show the mass splitting at all of our six quark 
masses. The reason that we are able to recover a mass splitting for an additional two lighter 
quark masses is because, as discussed in the introductory section, correlated errors in the 
correlation functions are suppressed in constructing the effective mass splitting. The state 
extracted with the pentaquark operator indicates the possibility of binding for the four 
smallest quark masses shown. 

In the infinite- volume limit, the lowest energy two-particle state, the P-wave N + K, 
approaches the energy of the .S-wave N + K state. As a genuine single-particle state is 
expected to have a small volume dependence on our lattice, we also show in Fig. ^] the 
mass splitting with the N + K two-particle threshold. In this figure, the mass difference is 
also negative, as is necessary for the presence of a bound state, suggesting that the presence 
of binding may prevail on larger lattices. Finally, we emphasise the possibility that the mass 
splitting may decrease as the light quark mass regime is approached, allowing the transition 
to a resonance at physical quark masses. High statistics studies at lighter quark masses 
would obviously be of considerable interest. 
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FIG. 14: Mass splitting between the I(J P ) = 0(| + ) state extracted with the NK* pentaquark 
interpolator and the mass of the P-wave N + K energy. 
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FIG. 15: Mass splitting between I(J P ) = 0(| + ) state extracted with the NK* pentaquark 
interpolator and the two-particle S"-wave N + K mass threshold. 

C. Negative parity isovector channel 

For completeness we also include an analysis of the isovector channel. First we present 
the results for the isovector, negative parity channel. In Fig. EI] we see that the effects of 
the spin projection for the largest quark mass are small. This may be understood by the 
presence of the S"-wave A + K* and iV + K* two-particle states, which both the spin-| and 
spin-| projected correlation functions should couple to. We fit the effective masses extracted 
from the spin-| and spin-| projected correlation functions, shown in Figs.Eland^l at time 
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FIG. 16: As in Fig- El but for the isovector negative parity channel. 
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FIG. 17: Effective mass of the I(J P ) = 1(^ ) pentaquark obtained for the NK* interpolator, 
Xnk*- The data correspond to m n ~ 830 MeV (squares), 700 MeV (circles), and 530 MeV (trian- 
gles). 



slices 20 — 30 and 18 — 21, respectively. 

The masses are presented in Fig. EH along with the corresponding mass extracted from the 
NK interpolator and the relevant two-particle states. The mass extracted from the spin-| 
projected correlation function is in good agreement with the mass of the N + K* two-particle 
state. Although this correlation function must have a contribution from both the A + K 
and N + K* two-particle states, we are probably accessing an admixture of these states. A 
correlation matrix analysis is required to separate this admixture. The mass extracted from 
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FIG. 18: As in Fig. EH but for the I{J P ) = state. 

the spin-| projected correlation function is consistently more massive in this channel than 
the mass extracted with the NK interpolator. However, in Ref. \^ we used a correlation 
matrix to extract the mass from the NK interpolator. This process removes excited state 
contamination and renders the ground state mass consistently smaller. In addition, by 
considering the two-point functions calculated with the NK and NK* interpolators at the 
quark level, one would naively expect the NK* interpolator to couple much more strongly 
to the N + K* two-particle state than the NK interpolator. In fact, the weak coupling of 
the spin-| interpolators to the lowest lying spin-| states is reflected in the relatively large 
statistical uncertainties for the NK* interpolator results. 

In summary, as in the case of the negative parity isoscalar channel, there is no indication 
of the lattice resonance signature. 

D. Positive parity isovector channel 

Finally we complete our discussion with the positive parity isovector channel. In Fig. 1201 
we see that the spin projection suggests the presence of two distinct states but the errors 
overlap. The effective masses corresponding to these correlation functions are presented in 
Figs. ED and E21 We fit the effective masses calculated from the spin-| and spin-| projected 
correlation functions at time slices 19 — 22 and 15 — 20, respectively. Due to the poor signal, 
results for only the three largest quark masses are shown. The masses of these states are 
presented in Fig. 1221 along with the mass extracted with the NK interpolator and energies 
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FIG. 19: Masses of the J( J p ) = and 1(§~) states extracted with the NK* interpolating field 

as a function of m^. For comparison, we also show the mass of the I(J P ) = 1(^ ) state extracted 
with the NK pentaquark interpolator in a correlation- matrix analysis a], and the masses of the 
S*-wave N + K, N + K* and A + K two-particle states. Some of the points have been offset 
horizontally for clarity. 
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FIG. 20: As in Fig. [2J but for the isovector positive parity channel. 



of the relevant two-particle states. Neither state extracted with the NK* interpolator lies 
below the lowest energy scattering states, which is necessary for binding. Therefore the 
analogous lattice resonance signature is also absent in this channel. 
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FIG. 21: Effective mass of the I{J P ) = 1(^ + ) pentaquark obtained from the NK* interpolator. 
The data correspond to ~ 830 MeV (squares), 700 MeV (circles), and 530 MeV (triangles). 
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FIG. 22: As in Fig. ED but for the I(J P ) = 1(f ) state. 



IV. CONCLUSIONS 



We have completed a comprehensive analysis of the isospin and parity states of the spin 



| pentaquark. Following our previous work we search for the standard signature of a 
resonance in lattice QCD, where the presence of attraction renders the resonance mass lower 
than the sum of the free decay channel masses at quark masses near the physical regime. 
This standard lattice resonance signature has been observed for every conventional baryon 



resonance ever calculated on the lattice 
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In the case of a pentaquark resonance, the relative mass splitting is expected to vanish in 
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FIG. 23: Masses of the I(J P ) = and 1(| + ) states extracted with the NK* interpolating 

field as a function of m^. For comparison, we also show the mass of the I(J P ) = state 
extracted with the NK pentaquark interpolator and the masses of the P-wave N + K, N + K* 
and the S'-wave N* + K two-particle states. Some of the points have been offset horizontally for 
clarity. 



the heavy quark limit. Therefore the analogous lattice resonance signature for a pentaquark 
state will exhibit a negative mass splitting at intermediate quark masses, with a general 
trend towards zero as the heavy quark limit is approached. 

In our examination of spin-| pentaquark states we have discovered evidence of the stan- 
dard lattice resonance signature, in the spin-| positive-parity isoscalar channel. At interme- 
diate quark masses, the presence of attraction between the constituents of the pentaquark 
baryon is sufficient to render the mass of the pentaquark state lower than the N + K two- 
particle threshold. The mass splitting approaches zero as the light quark masses become very 
large, in accord with expectations. Moreover, Fig. ITol suggests that the resonance signature 
might prevail on larger lattices, and provides a hint that the mass splitting will decrease as 
the light quark mass regime is approached, allowing the transition to a resonance at physical 
quark masses. 

Future work must explore the volume dependence of the binding observed in the spin-| 
positive-parity isoscalar channel. Otherwise, one cannot completely rule out the possibility 
that the observed binding is a pure finite-volume effect, reflecting a non-trivial scattering 
phase shift [29] in the I{J P ) = 0(| + ) N K scattering channel at momentum p = 2tt/L. If 
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this is the case, the mass splitting will go to zero as the volume of the lattice is increased. On 
the other hand, a genuine single-particle state is expected to have a small volume dependence 
such that the negative mass splitting is preserved in the infinite volume limit. As Fig. ^] 
illustrates, the observed mass splitting with the N + K two-particle threshold is sufficient 
to maintain the lattice resonance signature in the event that the volume dependence of the 
spin- 1 positive-parity isoscalar state is indeed small. 

The I{J P ) = 0(| + ) channel is therefore an interesting pentaquark resonance candidate 
for further study. Chiral fermions such as the overlap fermion action 30] allowing access 
to the lightest quark masses 0, [3| should be brought to bear on this particular channel. 
High statistics studies will be vital in rendering a conclusive result. Alternative resonance 
s lgnat ures such as the veW dependence of the residue QQc, „_e uuder hybrid 

nn 

boundary conditions |9|, [10J should also be brought to bear on this most promising channel. 
Ultimately, it will be important to investigate the nature of this state in full QCD, where 
the dynamical generation of quark-antiquark pairs is accounted for in the construction of 
the gauge field ensemble. 
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